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Abstract: We present explicitly Poisson structures, for both time-dependent 
and time-independent Hamiltonians, of a dynamical system with three degrees 
of freedom introduced and studied by Calogero et al [2005]. For the time- 
independent case, new constant of motion includes all parameters of the system. 
This extends the result of Calogero et al [2009] for semi-symmetrical motion. 
We also discuss the case of three bodies two of which are not interacting with 
each other but are coupled with the interaction of third one. 



1 Introduction 

In references [l],[2] Calogero et al. suggested and studied the system of differ- 
ential equations 
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as a model for motions of three bodies where forces determine velocities (hence 
the name Aristotelian), rather than accelerations . They were considered to 
be a prototype for a large class of models exhibiting transitions from simple 
to complicated motions that can be explained as travel on Riemann surfaces 
(see also extensive list of references in [T] and the fortcoming article [5]). In 
this work, following references [4] and [5] , we shall present Poisson structures of 
an equivalent system of differential equations obtained after elimination of the 
linear one-body forces in Eq.(IT]). 

Differentiation of Eq.([T]) implies that they may be interpreted as Newtonian 
equations for an inverse-cubic force field. Historically, cubic forces were first 
considered, in the Newtonian context, by Jacobi [6j. Poincare, in his search of 
infinitely many periodic solutions of three body problem, replaced the inverse- 
square force with an inverse-cubic one . The relation of inverse-cubic forces 
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with integrable three body problem of Newtonian type with inverse-square po- 
tential were explained in an appendix of reference [2]. See also [8]- [12] for other 
treatments of inverse-cubic forces. 

The system in Eq.((T]) was communicated to us by Yavuz Nutku who also 
presented its time-dependent conserved Hamiltonian functions [13] (c.f. Eqs.(|3])- 
A transformation of variables (c.f. Eq.©), first appeared in reference [T], 
takes one of the Hamiltonian functions into time-independent form but the 
second one still remains time-dependent. Using techniques of reference [3], we 
shall present, in section (3), a formal Hamiltonian structure with this time- 
dependent Hamiltonian function. 

We shall indicate, in section (4), a potential function by which the dynamical 
equations can be cast into a gradient flow. A linear change of variables that 
includes time-independent Hamiltonian function as one of the new coordinates 
reduces the number of equations into two which are still in gradient form. We 
shall then be able to integrate for the second time-independent conserved Hamil- 
tonian. For the full-symmetrical and semi-symmetrical cases, characterized by 
a = b = c and a = b ^ c, respectively, the bi-Hamiltonian form of equations 
for Aristotelian model of three body motion will be presented in section (5). 
It becomes necessary to give a separate treatment of a subcase of the semi- 
symmetrical motion characterized by a = & ^ and c = 0. This correponds to 
a model of three body motion in which two of the bodies do not interact with 
each other, yet their motions are coupled due to their interaction with the third 
body. 

In section (6), we shall construct only the second time- independent Hamil- 
tonian function containing all coupling constant, which is enough to cast the 
system into bi-Hamiltonian form in three-dimensions. We shall first complete 
the discussion for non-interacting two body case by considering a ^ b ^ and 
c = 0. We shall then present the general form of Hamiltonian function for 
generic values of coupling constants. This generalizes a result of reference [2] 
where a time-independent conserved function was found by a different analysis 
for the semi-symmetrical case. Present results are contributions with a geomet- 
ric taste to the ongoing investigations [3] of the authors of references [I] and [5] 
where they have been presenting full analysis of the dynamics such as methods 
of integrating Eq.(|T]), analytic structure of solutions, as well as the geometry of 
the Riemann surfaces involved. 



2 Aristotelian Model of 3-Body Motion 

In [1], the motion of three point particles, for which the velocities, rather than 
accelarations, are determined by forces, was modelled by the equations 

i I 9n-2 , gn-l 



where gn are coupling constants, z„ are (complex) coordinates of particle po- 
sitions, rt = 1, 2, 3, > 0, and i = \/— T. Ovcrdots denote derivatives with 
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respect to (real) time variable t and Eq. ^ was referred to as physical model. It 
was remarked that this systems admits canonical Hamiltonian formalism with 
the Hamiltonian function 

TTf \ -c^S r • I Pn+l ~ Pri+2 T /„N 

H{z,p} = T,,^^i[-lUJZnPn + gn J (3) 

Zn+1 — Zn+2 

where the momenta p„ look like Lagrange multipliers. 

We shall consider the system of equations ©, and not its extension by 
the addition of momenta, with zi = x, Z2 — y, z^ — z and gi = a, (72 = 
b, g^ — c resulting in Eq. ^ . Nutku 13| indicated two time-dependent conserved 
Hamiltonian functions 

H^^\x,y,z) = e-"^\x + y + z) (4) 

H^^Hx,y,z) = le-^^^'{x^+y^ + z^)-{a + b + c)t (5) 

for Eq.(IT]). The choice of variables, first appeared in [1], 

{u,v,w,T)=e-''^\x,y,z,^^—) (6) 

factors out the linear forces and transforms the system in Eq.([T]) into the one 
with two-body interactions only 
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Here, prime denotes derivative with respect to complexified time t and, for this 
reason Eq.Q was referred to as auxiliary model [1 . This transformation makes 
the first Hamiltonian function H^^^^ independent of time, but nevertheless i/*^^^ 
is still time-dependent 

H'^^\u,v,w) = u + v + w (8) 
H^'^\u,v,w) = ^v"^ + v? -2{a + b + c)t. (9) 

These relations were used in [2] to obtain the general solution of the system ([7]). 
Level sets of i?*^^) are planes and of i?*^^) are spheres expanding with velocity 



2(a + b + c) / y 2{a + b ^ c)t + where is the initial radius of the sphere, 
represents the coordinates of center of mass of three bodies. We note that 

H'^^^ (u, V, w) — uv + uw + vw + (a + b + c)t 

is also a time-dependent conserved function for Eq.([7]), but it is not independent 
from and 

i/(3) = i((i7(l))2_^(2))_ 



3 



Moreover, there is another relation between the time-dependent Hamihonians 

(m - vf + {v- wf + {w- uf = 2i7(2) _ 2H^^^ + 4(a + + c)t 

which can be used to ehminate the time variable r. This seemed to be useful in 
obtaining second time-independent Hamiltonian function in somewhat different 
analysis of reference [2] . 

Apart from the semi-symmetrical and the full-symmetrical (or integrable) 
configurations, the two-body case was obtained by the restriction a — b — 
and c 7^ pn which is not much interesting in the present framework. However, 
in presenting bi-Hamiltonian structures, we shall treat separately a particular 
case of semi-symmetrical motion in which the interaction between two of the 
bodies is neglected. Note also that, by a rescaling of time r one of the coupling 
constants can be removed from the right hand side of equations ([7]). In this 
work, we shall consider the variables u — {u,v,w) in a real domain. In the 
sequel ?7 = U • V = li • V will denote the vector field associated with Eq.©. 



3 Time-dependent Poisson structures 

We first present the Hamiltonian structure of equations ([7]) with the time- 
dependent Hamiltonian function H^'^K This can be achieved by considering 
Poisson structures in space-time variables (t,u). In other words, we add the 
equation f = 1 to the system in Eq.(I7]) (if necessary, after performing a linear 
change of parametrization by time) . The fiow of suspended system in four di- 
mensions is then generated by the vector field dt + U and our task is to find 
a Hamiltonian representation of this. We shall use contravariant description of 
Hamiltonian formalism by Poisson structures. 

A Poisson structure on a manifold N is defined by a skew symmetric con- 
travariant bilinear form subjected to the Jacobi identity expressed as the van- 
ishing of Schouten bracket of Poisson tensor with itself [I1]-[I7]. Following [3], 
for a Hamiltonian formalism on a time-extended space N — I x AI, / C M or C 
and M C ffi'^ for the present context, we take the bi- vector field 

A{t,u) = V{t,u) Adr+n{T,u) (10) 

where V and H are time-dependent vector and bi-vector fields on Af , respec- 
tively. The Jacobi identity for A can be computed in a coordinate independent 
way using identities of Schouten algebra of multi- vectors. 

Proposition 1 A is a Poisson bi-vector field on I x M if and only if 

dn dV 

[n,H]=2yA— , [U,V] = VA—. (11) 

The Hamiltonian form of the suspended vector field dt + U on I x M will 
then be 

dt + U = A{dH) = V{H)dt + n{dH) - HjV (12) 
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where H is a time-dependent conserved function of U. Due to non-linearity of 
Jacobi identity, finding bi- vector A from Eas. pA|) and ([T2|) for given Hamiltonian 
function H, is a difficult task. The next result, which is the linearization of 
Jacobi identity by Hamiltonian vector fields, may help. 

Proposition 2 For the Hamiltonian vector field in Ea. U2\) . the Poisson bi- 
vector \10\) satisfies the infinitesimal invariance conditions 

dV dli dU 

-^ + [u,v] = Q, ^ + [c/,n] = yA^ (13) 

which are equivalent to the Jacobi identity in Eg. Moreover, V is an in- 

finitesimal automorphism of A. 

This follows from the Jacobi identity (of Schouten algebra) for multi- vectors 
(A,A, /i). Conversely, one obtains Eq. pTj) by inserting Hamilton's equations 
(fT2|) into Eq. (IT3l) . The last conclusion follows from second of Eq. dTTI) . 

Thus, the construction of Hamiltonian structure of a system admitting time- 
dependent Hamiltonian function amounts to solving the linear system consisting 
of Eqs. dT^ . (|13p and the conservation law for H. As it can be inferred from 
Eq. p^ the function H is, in addition, coupled to the vector field V by the 
condition V{H) — 1. 

The construction of Hamiltonian structure for the Aristotelian model of dy- 
namical equations relies on the observation that the first of Eq. ([T3l) gives a 
characterization of the vector field as a time-dependent infinitesimal symme- 
try of U. Denoting 

E = udu + vdy + wdw 

the Euler vector field and, noting that U is homogeneous of degree —1, it is easy 
to see that the vector field 

^ rr , 2tc 2t6 , ^ , 2t a 2tc , 2r& 2Ta 

E~2tU = {u )du + iv )d^u + iw )du, 

U — V u — W V — W V — u w — u w — V 

is a time-dependent infinitesimal symmetry of U, that is, 

[dr + U,E -2tU] =0. 

This condition is essential to show that the bi- vector field 

A={E~ 2tU) Adr + EAU (14) 

is Poisson [A, A] — 0. We can then cast the system ([7]) with time-dependent 
conserved function in Eq.(|9]) into an autonomous Hamiltonian system in four- 
dimensions. 

Proposition 3 The vector field (1,U) is Hamiltonian for the Poisson bi-vector 
in Eq. and with the Hamiltonian function defined by the time- dependent 
function in Eq. ^ 



d^ + U = A{dH), H=-ln.H^^'> 
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This rather formal Hamiltonian structure may be useful in investigation of 
geometric structure of solution space as well as symmetries and invariants of the 
flow thereon [3]. For example, it follows that the time-dependent infinitesimal 
symmetry E—2tU of C/ is a Hamiltonian vector field for A with the Hamiltonian 
function r. In other words, the solution space of Eq.© may be realized as level 
sets of Hamiltonian function of an infinitesimal symmetry of motion it describes. 
The conserved function H^-^'' gives the Hamiltonian vector field H'^'^^dr + U) 
whereas results in —2H^^^U. 



4 Potential Function 

The force field of Aristotelian model can be derivable from a logarithmic poten- 
tial field. 

Proposition 4 The dynamical vector field U in Eq. ^ generating the Aris- 
totelian motion of three bodies is a gradient field U ='S/ F, with the potential 
function 

F = ln{v - w)''{u - wfiu - vf. 

Note that the potential function is a solution of the partial differential 
equation 

duF + dyF + droF ^ 

which is another manifestation of the conservation law for H'-^\ 

We shall consider a linear transformation of Cartesian coordinates (w, u, w) 
by which the conserved Hamiltonian function H^^^^ is eliminated and the po- 
tential function F becomes a function of two variables. The reduced system in 
two variables, which is still a gradient fiow, admits a non-canonical symplectic 
structure similar to the one for the Lotka-Volterra system in [18] . 

Proposition 5 In the orthonormal coordinates 

Q= ^{u + v + w), T] = ^{u~v), £_ = ^{u + v~2w) (15) 

the potential function F becomes a function of if], ^) only 

^^(ry,0=lnrf(V3e + r7)'(V3e-?7r- (16) 
In the new coordinates the right handed orthonormal basis vectors are 

ei = VC, 62 Vr;, eg = V^. 
From the inversion of Eq.(|15|) we observe that the combinations 



u — V — V^rj, V — w — —^{-\/3^ ~ v): w — u — —{Vs^ + 77) 

V 2 V 2 



are independent oi ~ H^^^^ 
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Proposition 6 On level sets of Hamiltonian function H^^^ with coordinates 
(77,^), the dynamical system in Eq.^ becomes 



c 



which is also a gradient system with the potential in Eq. il6]) . 

These are the reduced equations of motion after elimination of the motion 
of center of mass. General theory of reduction implies that the reduced space 
must be symplectic [T7]. In fact, as any orientable two dimensional space is 
symplectic, we can take an area form (f>{r],^)d^ A drj as the symplectic two- 
form for Ea. (IT7|) . In the context of local structure of Poisson manifolds, as 
described in reference [15], this reduction gives the symplectic foliation of the 
three dimensional space of variables (u, v, w). Liouville theorem for the adapted 
symplectic form and Eq. (|17p implies the invariance of the volume elements in 
the sense that the volume density (/'(t],^) satisfies 

where V denotes (9^, d^). The characteristics of this equation is the second time- 
independent conserved Hamiltonian we are looking for. It is also the Hamilto- 
nian function for Eq.(fT7|) with respect to the symplectic structure introduced 
above. 



5 Bi-Hamiltonian structures in 

The construction of bi-Hamiltonian structure in three dimensions requires two 
(time- independent) Hamiltonian functions Hi, H2 and a conformal factor (this 
will be seen to be the same function involving symplectic structure mentioned 
above) such that Hamilton's equations take the form 

u = (t>VHi X VH2. 

The Poisson tensors can be identified with vectors <j>WHi and — 0Vi?2 using the 
isomorhism between three vectors in M.^ and 3x3 skew-symmetric matrices. 
The corresponding Hamiltonian functions are H2 and Hi , respectively. Poisson 
tensors constructed this way can always be made into a compatible pair. A 
priori unspecified function (f) is related to an invariance property of the Jacobi 
identity in three dimensions. Namely, any multiple of a Poisson tensor with 
an arbitrary function is also a Poisson tensor. Thus, to cast a system into bi- 
Hamiltonian form, we need to find fundamental conserved quantities and then 
determine the multiplicative function <j>. See references [S],[Tn]-[17] for more 
details and various explicit examples. 

Since we have the time-independent conserved Hamiltonian i/*^^-' given by 
Eq.®, our first task is to find a second one. We want to search for the second 
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time- independent conserved quantity for the system in Eq.© or, equivalently, 
in Eq. (fT7|) using the reduced gradient flow of the latter. Note that, if there exist 
a time- independent conserved Hamiltonian function 7J(?7,5), then we can write 

VF = 0VC X Vi? = ei X (f)VH 

for some function </). It then fohows that the non-zero part of Ea. ([T7|) admit 
(noncanonical) symplectic formulation in the (77, ^)-variables 

d,^F = -<t>d^H, d^F ^ (j)dr^H (18) 

with the symplectic two- form 0(?7,^)dC ^ drj which is always closed and non- 
degenerate in two dimensions. Thus, the linear change of coordinates in Ea. (|15l) 
enables us to realize the local structure of Poisson manifold that we are going 
to construct. More precisely, the symplectic foliation of the space of variables 
(C, ry, ^) consists of coordinate planes C =constant, or equivalently, the level sets 
of the conserved Hamiltonian function H'^^\ 

The Hamiltonian function (77, ^) will be a function of characteristic curves 
defined by the Hamiltonian system in Eq. ljlSp and it may depend arbitrarily on 
the variable C- To find the characteristic curves of Ea. (|T8| . or what we shall 
call the fundamental conserved quantity, we eliminate the time derivatives in 
Eq. pTp and obtain the ordinary differential equation 

dri {a + b + c)r]^ + V3{a - b)rj^ - 3cf 

which is homogeneous of degree two. So, for 6 = ^/-q, its solution defines the 
characteristics 

constant = In 77+ / i— ^ ^ ^.dO (20 

J V3(a -b) + (4a + 46 + c)e + \/3(a - 6)6*^ - icO^ 

or, equivalently, the fundamental conserved quantity for the flow. That means, 
the time-independent Hamiltonian function we are seeking is a function of C and 
the function defined by Ea. (l20t . Starting from the simplest, we shall first present 
explicitly bi-Hamiltonian structures of symmetrical cases and then proceed, in 
the next section, to analyse the generic solution for the fundamental conserved 
quantity in Ea.(P(7|. We shall see that, in semi-symmetrical motion, the subcase 
c = requires a separate treatment and we shall do this in the last subsection. 

5.1 Pull-symmetrical case 

We set a = 6 = c for which Eq. (fT9| becomes independent of coupling constants 

d^ _ 277C 
dry if - ^2 

and the solution gives 



(21) 



hf{(,,v)=e-Uv' (22) 
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which can easily be verified to be conserved under the flow of Eq. ([T7| . Simple 
manipulations show that 

0/(^,7?) = l/V3ry(3f-772) 

is the conformal factor entering the definition of symplectic structure. As a 
function of variables {u,v,w) in Eq.Q, we find 

Proposition 7 For the Aristotelian model with equal coupling constants, the 
dynamical equations 
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admit bi-Hamiltonian structure 

u = P/i(u)Vi7/(u) = P/2(u)Vi?(i)(u) 
with the following pairs of Poisson matrices and Hamiltonian functions 

1 / 1 -1 

V6[u - v)[v - w)[w ~ u) y 1 -1 
Hf{u) = -^{u + v -2w)[{u + v -2wY -9{u-v)'''], 
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The conformal factor for these Poisson structures is the function i^j(u) given 
as the coefRcient of the constant matrix of P/i(u). This function constitutes an 
invariant volume density for the solution space of variables u in the sense that 
the three- form (j)j{u)du A dv A dw has vanishing Lie derivative with respect to 
the vector field defined by the right hand side of Eq. (| 

5.2 Semi-symmetrical case 

Setting a ~ b ^ c, we have 

d^ 6a?7^ 



dr] (2a + c)t]^ ~ 3cC 



(24) 
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and its integration gives 



vi-Vii-f - ^^l*'"'''^' = constant (25) 

where, following reference [5], we introduce the constant ^ = {2a + c)/{8a + c). 
Various powers of the function in Eq. ([25)) can be adapted as the fundamental 
conserved quantity. In any case, we have to exclude the values fi — 1/A (or 
c = 0) and /i = 1 (or a = 0) in the following discussions. For a = 0, the 
evolution in the variable w disappears (c.f. Eq. ()26p ) and we are left with a two 
dimensional system. For c = 0, we have the special case of three body motion 
in which two of the bodies do not interact with each other. This necessarily 
requires separate treatment which we will take up in the next subsection. We 
refer to extensive discussion in reference [2] where it was shown that values of 
the constant n, in particular, the real rational values, play important role in 
determining dynamical evolution of the model. 

For the purposes of having a generalization of the full-symmetrical case and 
the same invariant volume density, we choose the fundamental conserved quan- 
tity for semi-symmetrical case to be 

hs{^,v) = e'-^^'-'He - A* ^ 1,1/4 

so that, when a = c, or equivalently, fi — 1/3 we obtain the function ft./(C, v) in 
Eq.' 



Proposition 8 For the Aristotelian model with two coupling constants a 
and c 0, the dynamical equations 



u — V u — w 

a c , , 
-f 26 

V — W V — u 
a a 



w = 1 

w — U W — V 

admit Hamiltonian structure with the Poisson tensor 



" 2 {u — v){v — w){w — u) 



and with the Hamiltonian function 



1 -1 
-1 1 

1 -1 



H,{u) ^{u + v- 2u;)2''/(i-^) [{u + v- 2wf - ^ (u - v) 
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The second Hamiltonian structure is defined by the Hamiltonian function H^^^ 
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u + v + 1^ and the Poisson tensor 

c 1 ,v — w w — u 
12u — V w — u 



c 3li u — v 

+ ^ 



+ 



12 4/i - 1 (w - w){w 
c 1 — fl . 1 




Note that, in the limit /x=l/3( or a = c) both ^^(u) (defined as the 
multiphcative factor in -Pi(u)) and -ffs(u) reduce to some constant multiples of 
<pf{u) and Hf{u), respectively. 

5.3 Non-interacting two body case 

We consider a = b ^ and c = 0, that is, /x = 1/4. This case corresponds to 
a situation where two of the three bodies do not interact with each other but 
with the third one only. Or, to a situation where the distance |u — t;| is so large 
that the interaction between bodies at u and v can be neglected. Setting the 
only constant a = 1, we have 

Proposition 9 Aristotelian equations of motion 

1,1,1 1 

U = , V = , w = 1 

U — W V — w w — u w — v 

for two bodies at positions u and v interacting with a third one at w are bi- 
Hamiltonian with 

fnl(u) 
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iJ„(u) = ^{u + V - 2w){u - vf 
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6 {v — w){'w — u) 
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2 w — u V — w 

H^'^\n) = u + v + w. 

There remains the case of non-interacting two bodies each of which interact 
with the third one with different coupling constants a, 6, a ^ b. We will give 
the form of Hamiltonian function of this general non-interacting two body case 
in the next section. 
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6 General Form of Hamiltonian Function 



We have seen that casting the dynamical system modelUng Aristotehan motion 
of three bodies into bi-Hamiltonian form requires first the integration of Eq. (|20p 
to find the time-independent Hamiltonian function and then, the invariant vol- 
ume density. In this section, we shall discuss the general form of characteristics 
and their domain of validity. Our discussion will, by no means, be exhaustive. 

We first give the form of Hamiltonian function for motions of three bodies 
two of which are not interacting with each other. The remaining third body 
interacts with them with different coupling constants a and b. Integration in 
Eq.daOl) gives 



Proposition 10 Let two non-interacting bodies at positions u and v interact 
with the third one at w. Then, the Aristotelian equations of motion 



b , a , b 
u = , V = , w — 



u ~ w V — w w — u w — V 

admit the conserved Hamiltonians H^^'^ and 

u ~ V 



/ 

H{u) = 2V4fc2-lln- 



V2 

+ ( V4fc2 - 1 - fc) ln( " ^ ~ - V4fc^ - 1 + 2fc) (27) 
V3('u — v) 

u + V — 2w 



-^(v/4fc2 - 1 + k) ln(^ + v/4fc2 - 1 + 2k) 

V3(w — v) 

where we introduce the constant k — {a + b)/{y/i{a — h)). 

Returning to the general case, we first cast Eq. ipn]) into the form 

and, after some manipulations, obtain 

Proposition 11 Let 9± be the roots of quadratic term in numerator and, let 
0i, i — 1,2,3 be the roots of cubic polynomial in denominator of the integral in 
Eg. \20\]. Then, the second fundamental (time-independent) conserved quantity 
for Aristotelian model of three body motion is 

u — V 





93){di-( 
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To this end, we want to discuss examples of restrictions on the domain of 
definition of the fundamental Hamiltonian function h. The integral in Ea. (|28l) 
is singular for a real root ^rcai of the cubic polynomial in denominator. That is, 
for the case 9 = ^/r/ = ^loai- In the real variables {u,v,w) this implies 

U{1 - VSflroal) + V{1 + Vaflrcal) ~ 2w ^ 0. 

The cases described by the conditions u = v (or rj — 0), v — w and w = u are 
particularly included in this equation. Moreover, for non-zero values of ^rcai, 
by adding and subtracting u and v to the above equation we can obtain the 
singular cases 

V — W = 0, 1 — \/30rcal = 
W — U = 0, 1 + \/30roal = 

which let the variables u and respectively, be free but bring restrictions on the 
coupling constants a, 6, c. Finally, the present form of roots of cubic equation 
imposes some conditions on the coupling constants related to the discriminant 
of cubic polynomial. 

Proposition 12 For the values of the vector field u not perpendicular to the 
constant vector 

(l-\/30,.oal,l + \/30rcal,-2), 

the dynamical system in Eq.^ admits hi- Hamiltonian structure with Hamilto- 
nian functions in Eqs. ^ and I129\) . 

In the particular value ^?icai = of the real root, that may occur for example 
in semi- and full-symmetrical configurations, the condition of above proposition 
simply prevents the adapted coordinate ^ to become zero. For generic values of 
^rcai it just requires us to be away from the line ^ = OicaiV on level sets of H^^\ 

The discussion on domain of validity of bi-Hamiltonian structure will become 
more conclusive if we adapt the parameters 

p a-b a+b 

for the integration of Eq. (pn| . In these parameters, the semi-symmetrical case 
is characterized by 

2a I 1- ^^ 

p — 0, a = — = 

^ ' ^ 3c 3 - 1 

and for the full-symmetrical case we have p — 0, q — 2/3. Note that, the 
non-interactive two body case (c = 0) is necessarily excluded from the present 
discussion. The integral for the fundamental conserved quantity takes the form 

/iff, 77) = In 77+ / — ^-^-^ :rd-d (30) 

^ ' 7 -2p(p2 + i8(7+15)/27-((p2 + 12g+l)/3)7? + 7?3 ^ ^ 
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which can be put into the form of Eq. (pS)) with the roots 

^1 = ---^(l+p2 + 12g + A2/3) 

^± = --^(^±(l+p' + 12g) + ^^A2/3) 
GA^'"' 

of the cubic denominator. Here, ^± = 1 ± i-\/3 and 

A = -pip^ + 18q + 15) + ^/27p^ + 6p2(6(13 - 3q)q + 37) - (1 + 12q)\ 

In the semi-symmetrical case, we have A = (—1 — I2q)^^^, q = 2a/3c and so 
di = 0, 'd± — y/1 + 12q/ a/S.Iu the full-symmetrical case, q — 2/3, A = 27i and 
■d^ = 0, i9± Vs. 

As a final example of this analysis, we want to present the connection be- 
tween values of the discriminant and the constants appearing in the integral of 
Eq. (pO| . Recall the discriminant 

A = 4(-27/ + 6p2(_37 -I- 6g(-13 + 3q)) + (1 + 12g)3)/27 (31) 

of cubic polynomial in Eq. pop which is the negative of the square- rooted term 
in A. Regarding the right hand side as a quadratic polynomial in p^ , we find its 
discriminant to be 

- 16777216(1 + 3(?)2(7 + 3qf{l + 12^)^177147. (32) 

A will be a perfect square at the values oi q = —1/3,— 3/7,— 1/12 and these 
correspond to the following sets of values of the coupling constants 

a + b + c^O, 7{a + b) + 9c^0, 4(a + 6) + c = 

two of which are seen in Eq. (j20l) . First two restrictions imply, for the semi- 
symmetrical motion, fi — and fj, = 1/29, respectively, the last one, however, is 
not possible. 

We expect the results of the present work to be useful in an exhaustive 
analysis of various cases of Aristotelian model of three body motion as well as 
in understanding the physics behind this model. 
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